Charged-current weak interaction processes in hot and dense matter and its impact 
on the spectra of neutrinos emitted from proto-neutron star cooling 
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We have performed three-flavor Boltzmann neutrino transport radiation hydrodynamics simula- 
tions covering a period of 3 s after the formation of a protoneutron star in a core-collapse supernova 
explosion. Our results show that a treatment of charged-current neutrino interactions in hot and 
dense matter as suggested by Roddy et al. [Phys. Rev. D 58, 013009 (1998)] has a strong impact 
on the luminosities and spectra of the emitted neutrinos. When compared with simulations that 
neglect mean field effects on the neutrino opacities, we find that the luminosities of all neutrino 
flavors are reduced while the spectral differences between electron neutrino and antineutrino are 
increased. Their magnitude depends on the equation of state and in particular on the symmetry en- 
ergy at sub-nuclear densities. These modifications reduce the proton-to-nucleon ratio of the outflow, 
increasing slightly their entropy. They are expected to have a substantial impact on the nucleosyn- 
thesis in neutrino-driven winds, even though they do not result in conditions that favor an r-process. 
Contrarily to previous findings, our simulations show that the spectra of electron neutrinos remain 
substantially different from those of other (anti)neutrino flavors during the entire deleptonization 
phase of the protoneutron star. The obtained luminosity and spectral changes are also expected to 
have important consequences for neutrino flavor oscillations and neutrino detection on Earth. 
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Supernova explosions of massive stars are related to 
the birth of neutron stars due to the collapse of the stel- 
lar core at the end of stellar evolution [2]. The detec- 
tion of neutrinos from SN1987A [3] confirmed that the 
ps 3 x 10 53 ergs of gravitational energy gained by the core 
collapse are emitted as neutrino radiation on time scales 
of tens of seconds, during which the central protoneutron 
star (PNS) cools, dclcptonizes and contracts. In the de- 
layed neutrino-heating explosion mechanism [2, 4], neu- 
trinos also transport energy from the nascent PNS to the 
stalled bounce shock. This mechanism remains the most 
viable scenario to explain supernova explosions as con- 
firmed by recent two-dimensional simulations [5]. Once 
the explosion sets in, the continuous emission of neutri- 
nos from the PNS drives a low-mass outflow known as 
neutrino-driven wind [G] that is currently considered the 
favored site for the productions of elements heavier than 
iron (e.g. [7]). As neutrinos travel through the stellar 
mantle, they can suffer flavor oscillations [8], contribute 
to the nucleosynthesis of several rare isotopes [9] and even 
drive an r process in the He-shell of metal-poor stars [10] 
before they are finally detected on Earth. 

Accounting for all aspects discussed above requires the 
knowledge of the spectra of the neutrinos emitted during 
the cooling phase of the PNS. Due to their low ener- 
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interact only via neutral-current reactions. 



Hence, together with the neutron-richness of the PNS 
surface one expects the following neutrino-energy hier- 
archy: (E^ T ) > (E„ c ) > {E Ve ) [11, 12], with (E) the 
average energy of the neutrino spectrum. Early super- 
nova models [13] predicted large energy differences be- 



tween v e and v e that resulted in neutron-rich ejecta as 
required by r-process nucleosynthesis [14]. However, as 
the treatment of neutrino transport and weak interaction 
processes improved, the computed energy difference be- 
tween v e and v e decreased and the early wind ejecta be- 
came proton rich [15]. More recently, it has been possible 
to perform supernova simulations based on three-flavor 
Boltzmann neutrino transport for time scales of several 
tens of seconds [16, 17], covering the whole deleptoniza- 
tion of the PNS. These simulations predict a continuous 
decrease of the energy difference between neutrinos and 
antincutrinos of all flavors that became practically indis- 
tinguishable after ps 10 s. The exact value depends on the 
progenitor. The proton-richness of the ejecta increases 
continuously with time and leaves the vp process [18] as 
the only mechanism for producing elements heavier than 
iron in neutrino-driven winds. 

The simulations of ref. [17] have been recently ana- 
lyzed, showing that the convergence of neutrino and an- 
tineutrino spectra at late times is due to the suppression 
of charged-current processes at high densities [12]. This 
analysis was based on a set of neutrino opacities that as- 
sume a non-interacting gas of nucleons and nuclei. This 
approximation may be valid during the accretion phase 
prior to the onset of the supernova explosion when the re- 
gion from where neutrinos decouple, the ncutrinosphcrcs, 
is located at relatively low densities, ~ 10 11 g cm~ 3 . 
However, as the PNS cools the neutriospheres move to 
increasingly higher densities where the non-interacting 
gas approximation breaks down. The nuclear interaction 
is treated at the mean-field level in equations of state 
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(EoS) commonly used in core-collapse supernova simula- 
tions [19, 20]. However, its influence on weak interaction 
processes is often neglected. In this Letter, we show that 
a treatment of neutrino-matter interactions that is con- 
sistent with the underlying EoS has a strong impact on 
the spectra and luminosities of the emitted neutrinos. We 
discuss the relevance for nucleosynthesis, neutrino oscil- 
lation studies and neutrino detection. 

Our study is based on the work of Ref. [1] where cor- 
rections to the opacities due to strong interactions are 
considered at the mean-field level. Effects of many-body 
correlations [21] will be considered in a forthcoming pub- 
lication. They are expected to affect the neutrino spectra 
at later times [22] than considered in the present study. 
We focus on charged-current (anti)neutrino absorption 
processes on neutrons and protons and the inverse reac- 
tions: e~ + p n + v e and e + + n p + D e , which are 
those mainly affected by mean-field corrections. 

EoS commonly used in core-collapse supernova simula- 
tions, see e.g. Refs. [19, 20], treat neutrons and protons as 
a gas of quasi-particles that move in a mean-field single- 
particle potential U. Assuming non-relativistic nucleons, 
which is a good approximation for densities p < 5po 
where po « 2.5 x 10 14 g cm -3 , the energy momentum 
relation closely resembles the non- interacting case [1]: 
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with particle rest- masses to,;. Both the single-particle 
mean-field potentials and the (Landau) effective masses, 
to* depend on density, temperature and proton-to- 
nucleon ratio, Y e . Importantly, due to the extreme 
neutron-rich conditions the mean-field potentials for neu- 
tron and protons can be very different with their relative 
difference U n — U p directly related to the nuclear symme- 
try energy [1]. 

In order to quantify the impact of the mean field po- 
tentials, let us consider (anti)neutrino absorption on neu- 
trons (protons). Assuming zero momentum transfer, i.e. 
Pn ~ Pp (clastic approximation), the electron(positron) 
and (anti)ncutrino energies are related by: 

E Ve = E e - - (m n - m p ) - (U n - U p ), (2a) 
= E e + + (m n - m p ) + (U n - U p ). (2b) 

Eqs. (2a) and (2b) show that the contribution of the 
mean field potential reduces (increases) the energy of 
the emitted neutrinos (antineutrinos) . The energy dif- 
ference between neutrinos and antineutrinos is increased 
by an amount 2{U n — U p ). The opacity, or inverse mean- 
free path, for (anti)neutrino absorption also changes. In 
the elastic approximation and assuming extreme rela- 
tivistic electrons, the opacity for neutrino absorption be- 



comes [1, 23]: 
x(E„ e ) = 
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with E Vc and E e related by equation (2a). The emis- 
sivity, j(E Uc ), can be obtained from the detailed bal- 
ance condition j{E Uii ) = exp(—f3(E u<; — p1 q ))x(E^ ), with 
— Me ~ (Mn ~ Mp) the equilibrium neutrino chemical 
potential, p the chemical potential including rest mass 
and /3 the inverse temperature. The opacity and emissiv- 
ity for antineutrino absorption are obtained exchanging 
neutron and proton and using equation (2b) to relate the 
positron and antineutrino energies. In equation (3), Gp 
is the Fermi coupling constant, V u d is the up-down en- 
try of the Cabibbo-Kobayashi-Maskawa matrix, gy and 
qa are vector and axial-vector coupling constants and 
n Pi „ the number density of protons or neutrons. / is the 
Fermi-Dirac distribution function and rj is the chemical 
potential (without rest mass). The quantity ip = r\ — U 
represents the chemical potential for a non-interacting 
gas of nucleons, that is related to the nucleon number 
density by the relation: 
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Eq. (4) provides a method of determining the mean-field 
potential, U, when using an EoS that does not provide 
this quantity, e.g., the EoS of ref. [20]. 




30 60 90 

Neutrino Energy (MeV) 



30 60 90 

Antineutrino Energy (MeV) 



FIG. 1. (Color online) Opacity and emissivity for neutrino 
(left panels) and antineutrino (right panels), evaluated at con- 
ditions p = 2.1 x 10 13 g cm" 3 , T = 7.4 MeV and Y e = 0.035. 



Figure 1 shows neutrino and antineutrino opacities and 
emissivities evaluated at conditions found at the antineu- 
trinosphere for the 18 Mq model of ref. [12] at 1 s after 

include the 
= -7.6 MeV 



bounce. The curves labeled RMF (U n ,U p 
contribution of the mean-field potentials U n 



and U D 



-14.2 MeV [20], while the curves labeled 
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(U n = Q,U P = 0) assume a non-interacting gas of nu- 
cleons, i.e. neglect the contribution of the potentials but 
still use chemical potentials as given by the EoS. Due 
to the presence of the mean-field potentials the effec- 
tive Q-value for electron capture increases with respect 
to the free case producing neutrinos with substantially 
lower energy. For the inverse process, neutrino absorp- 
tion, the opacity is enhanced due to the fact that the 
produced electron gains an energy U n — U p reducing the 
final-state Pauli blocking of the electron. The situa- 
tion is completely analogous to (anti) neutrino emission 
and absorption on heavy neutron- rich nuclei [24]. Using 
Eq. (3), it can be shown that the opacity for the non- 
interacting gas, XnU is related to the mean- field opacity, 
Xmf by Xmi{E) = xm{E + U n - Up). This relationship 
produces a large enhancement of the neutrino opacity 
at high densities, p ~ 10 14 , where U n — U p « 50 MeV 
when compared with the non-interacting approximation 
used in ref. [12]. For antineutrino absorption, due to 
the fact that positrons follow Boltzmann statistics, the 
non-interacting emissivity and mean-field emissivities are 
related by j m f(E) — j n ; (E — U n + U P ). The mean- field an- 
tineutrino opacity is larger at high densities as final-state 
Pauli blocking of the neutrons becomes less efficient. 

In the following, we explore the impact that a de- 
scription of opacities consistent with the EoS has on 
the spectra and luminosities of the emitted neutrinos. 
We have performed core-collapse supernova simulations 
based on spherically symmetric radiation hydrodynamics 
with three- flavor Boltzmann neutrino transport. Since 
our goal is to explore the differences in neutrino energies 
and luminosities due to the inclusion of mean- field poten- 
tials, we have used a low resolution transport scheme with 
12 energy bins and allowed only for radially in and outgo- 
ing neutrinos. Despite of its limited resolution, it repro- 
duces the absolute values of luminosities and average en- 
ergies predicted by higher resolution simulations [16, 17]. 
Table 1 of [12] list the weak processes considered in our 
simulations. We use the baryonic high-density EoS from 
Shen et al. [20] for matter in nuclear statistical equilib- 
rium (NSE) at temperatures above 0.45 MeV. As the 
tabulation of Shen et al. does not provide the mean-field 
potentials, we have computed them using eq. (4). In the 
non-NSE regime, we use the EoS of ref. [25], which we 
also use for electrons, positrons and photons in the NSE 
regime. The simulations arc based on the 15 Mq progen- 
itor of ref. [26]. Because spherically symmetric simula- 
tions do not result in explosions for such a massive iron- 
core progenitor, we enhance the neutrino heating rates 
in the gain region following the scheme of ref. [17]. It 
results in the onset of explosion at about 350 ms post 
bounce. The simulations are evolved from core collapse, 
through the explosion up to more than 3 seconds after 
bounce. During core collapse and post-bounce accretion 
phases, the mean-field potentials, U n , U p , are only on the 
order of several 100 keV in the region of neutrino decou- 



pling, which is located at intermediate densities on the 
order of 10 11 g cm~ 3 . Hence, the inclusion of mean-field 
potentials does not affect the supernova dynamics prior 
to the explosion. 
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FIG. 2. (Color online) The left panel shows the neutrino 
luminosity (upper) and average energy (lower) evolution. The 
right panel shows the evolution of the proton-to-nucleon ratio, 
Y e , and the entropy per nucleon for several mass elements 
ejected from the PNS surface. The curves shown in blue use 
neutrino opacities computed using the mean-field potentials 
of the EoS [20] while they are neglected on the red curves. 

After the onset of the explosion the neutrinospheres 
move to increasingly higher densities reaching values of 
the order of 10 13 g cm -3 . The left panels of Fig. 2 show 
the evolution of the luminosity and average neutrino en- 
ergy for all neutrino flavors. These observables are sam- 
pled in a co-moving reference frame at a distance of 
1000 km. Using charged-current neutrino opacities that 
include the mean-field potentials slightly reduces the lu- 
minosities for all neutrino flavors. Moreover, as expected 
from the discussion above, it enhances the differences in 
luminosities and average energies between neutrinos and 
antineutrinos. 

Fig. 3 shows the different neutrino spectra for all fla- 
vors at 3 s after bounce at a distance of 30 km outside the 
neutrinospheres. At this distance neutrinos can be con- 
sidered free streaming but they have not yet been subject 
to collective neutrino flavor oscillations [8]. These may 
result in spectral swaps [27] that occur in regions near to 
spectral crossings. We expect substantially different os- 
cillation patterns for the spectra obtained with opacities 
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FIG. 3. (Color online) Neutrino spectra for all flavors (solid 
lines left panel: u e , dashed lines left panel: ^ M / T , solid lines 
right panel: dashed lines left panel: P M / r )> including the 
mean- field potentials (blue) and without (red). 



considered in the present study. Furthermore, it is impor- 
tant to explore the sensitivity of the results to different 
EoS and in particular to EoS that are consistent with 
recent constrains on the nuclear symmetry energy [30]. 
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consistent with the EoS. 

The changes in electron (anti)ncutrino spectra and lu- 
minosities have important consequences for nucleosyn- 
thesis in neutrino-driven winds. The increased difference 
between average energies of v e and v e spectra impacts 
the Y e of the ejected matter resulting in slightly neutron- 
rich conditions for the early ejecta (see lower-right panel 
of Fig. 2) while at later times the ejecta become proton 
rich. In the simulation that neglects the contributions of 
the mean-field potentials the ejecta are always proton- 
rich. The decrease in Y e is not large enough to favor an 
r-process but may help in the production of isotopes like 
92 Mo that are only made under slightly neutron-rich con- 
ditions [18, 28]. Also relevant for nucleosynthesis is the 
slight increase in entropy per nucleon of the ejected ma- 
terial (see upper-right panel Fig. 2) that can be related 
to the reduced neutrino luminosities [2!)]. 

We have shown that a treatment of the charged-current 
(anti)neutrino opacities, that is consistent with the EoS 
as suggested by [1], has important consequences for the 
neutrino-spectra evolution during the PNS cooling phase. 
The most relevant finding is an increased difference be- 
tween average energies of v e and v e that persist during 
the whole simulation time of 3 seconds after the onset of 
the explosion. The changes on neutrino spectra are ex- 
pected to have important consequences for nucleosynthe- 
sis, flavor oscillations and neutrino detection on Earth. 
Our results imply that not only the evolution of the neu- 
trino luminosities [22] but also the spectral differences 
between v e and v e are sensitive to the symmetry energy 
of nuclear matter. Our simulations are based on neutrino 
opacities computed using the elastic approximation that 
neglects momentum exchange between nucleons. They 
need to be extended to consider the full kinematics [1] 
of the reaction and many-body correlations [21] that are 
expected to become important at later times than those 
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